Fractional topological phase for entangled qudits 
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We investigate the topological structure of entangled qudits under unitary local operations. Dif- 
ferent sectors are identified in the evolution, and their geometrical and topological aspects are 
analyzed. The geometric phase is explicitly calculated in terms of the concurrence. As a main 
result, we predict a fractional topological phase for cyclic evolutions in the multiply connected space 
of maximally entangled states. 
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In a seminal work, M. Berry 1] showed the impor- 
tant role played by geometric phases in quantum the- 
ory. Since then, the interest for geometric phases was 
renewed by potential applications to quantum computa- 
tion. The experimental demonstration of a conditional 
phase gate was provided both in Nuclear Magnetic Reso- 
nance (NMR) [H and trapped ions Q . Optical geometric 
phases have already been discussed both for polarization 
0] and vortex mode transformations [H, Q . The role of 
entanglement in the phase evolution of qubits was inves- 
tigated in refs.0,[l]- Recently, P. Milman and R. Mosseri 
[9l[l0| investigated the geometric phase and the topolog- 
ical structure associated with cyclic evolutions of arbi- 
trary two-qubit pure states. This structure has been ex- 
perimentally evidenced in the context of spin-orbit mode 
transformations of a laser beam 11 1 and in NMR [l2j ]. 
Although the topological nature of the phase acquired by 
maximally entangled states is well settled, the distinction 
between geometrical and topological phases has not been 
established clearly for partially entangled states. In this 
work we present a group theoretical approach which al- 
lows for a clear distinction between the two aspects. As a 
bonus, this approach is easily extended to higher dimen- 
sions, bringing an interesting prediction of a fractional 
topological phase. 

Let \ip) = - —1 0£n\ij) be the most general two- 
qudit pure state. We shall represent this state by the 
d x d matrix a whose elements are the coefficients ctij. 
With this notation the norm of the state vector becomes 
(ip\ip) = Tr(a' f a) = 1 and the scalar product between 
two states is (<p\ip) = Tr(/3'a), where (3 is the d x d 
matrix containing the coefficients of state \<j)) in the cho- 
sen basis. We are interested in the phase evolution of 
the state \xf>) under local unitary operations. So let us 
take two unitary matrices Ua and Ub belonging to U{d) 
and representing the operations performed in each sub- 
system separately. Under these unitary operations the 
state matrix will evolve as a(t) = Ua cc(0) U b , where 
Uj{t) = e^^Ujit) (j = A, B) and U 3 € SU{d). One 
can identify the following invariants under local unitary 
evolutions: Tr[pj], p = X,...,d, where pj is the reduced 
density matrix with respect to qudit j {pa = a 1 a* and 
Pb = eta'). In fact, the invariants are j-independent. 



The first one (jp = 1) is simply the norm of the state 
vector. One can readily relate the second invariant to 
the I-concurrence of a two-qudit pure quantum state 13 1 
C = \/2(l — Tr[p 2 ]), so that its invariance expresses the 
well known fact that entanglement is not affected by local 
unitary operations. The p = d invariant can be rewritten 
in terms of the former and T> = | det[a]|. In particular, 
for qubits we have C = 2 V. 

In the case of a cyclic evolution, UA(T)a(0)U B (r) = 
e lA ^a(0) . By taking the determinant of both sides we 



get: e la ^i> = 1 as long as V ^ 0. This implies that the 
possible acquired phases due to the SU (d) part of a cyclic 
evolution are A</> = 2nn/d, with n — 0, 1, 2, d — 1. 
For qubits (d — 2) one recovers the well known result 
A<p = 0, 7r. However, for d > 2 one obtains fractional 
phase values in steps of 2n/d. Now, we are interested 
in discussing in what sense this fractional phase can be 
considered as topological. For this aim, we will analyze 
the topology of the space of two-qudit states and how 
the total phase is built. In this regard, we would like to 



underline that according to ref . (14| , the geometric phase 
acquired by a time evolving quantum state a(t) is always 
defined as 



cf> g = arg<V(0)|V(*)> +* / dt M*)M*)> 



(1) 



that corresponds to the total phase minus the dynamical 
phase. Therefore, a topological phase, that is, an object 
that only depends on a given class of paths, can only find 
room as a part of the geometric phase, an object that is 
invariant under reparametrizations and gauge transfor- 
mations. Gauge invariance corresponds to the fact that 
the phase factors <fij(t) do not contribute to <f> g , which 
is completely determined by Uj(t), the sector where the 
fractional values occur. 

In order to characterize the space of states, we note 
that any invertible matrix admits a polar decomposition 
a = QS, where Q — \/V e M is a positive definite Her- 
mitian matrix, M is a traceless Hermitian matrix, and 
S = e**S, S £ SU(d). Since det[e M ] = e Tr ^ = 1, one 
easily finds det(a) = T)e ld ^. We can identify the time 
evolution as occurring in different sectors 



a(t) = <TDe 



i<t>{t) M(t) 



S(t), 



(2) 



2 



where we have denoted, cf)(t) — 0(0) + (pA(t) + 4>B{t) : 
M(t) = U A {t)M{Q)U A {t)\ and S{t) = U A {t)S{0)Ul(t). 
Therefore, we identify the evolution in three sectors of 
the matrix structure: an explicit phase evolution <p(t), 
an evolution closed in the space of traceless Hcrmitian 
matrices M(t), and the evolution S(i) closed in SU(d). 

Now we are able to discuss the topological aspects 
of the entangled state evolution in terms of these sec- 
tors. The space of positive definite Hcrmitian matrices 
Q has trivial topology. This is a noncompact manifold 
isomorphous to R d _1 , as it can be parametrized in the 



form Q 



AT. 



where f3 a are real numbers and T a 



(a = l,2,...,d 2 — 1) is a basis in the space of Hcrmi- 
tian traceless matrices. These T a 's are the generators of 
SU(d), so that 5 — e lu)aTa , with oo a real. They can be 
normalized in the form tr (T a T b ) = ^S ab and obey the 
Lie algebra [T a ,T b ] = if abc T c , where f abc are the struc- 
ture constants of SU(d). The first homotopy group of 
SU (d) is also trivial, however, the physical equivalence 
of a matrices differing by a global phase corresponds to 
considering the identification in SU(d), e l27rn/d S = S. 
This can be naturally implemented by associating the 
SU(d) sector of the matrix a with a corresponding sector 
for the quantum states, represented by transformations 
R(S) in the adjoint representation ST a S~ 1 — h a ■ T, 
h a = R(S)e a . In this manner, the matrices e i27Tn / d S are 
mapped to the same point R(S). In other words, a part 
of the evolution can be parametrized as R(t) G Adj(<i), or 
equivalently, in terms of a time dependent frame n a (t). 
Note that for qubits the adjoint representation corre- 
sponds to 5*0(3), the manifold used in ref. 0, [13] to 
describe maximally entangled states. An evolution S(t) 
starting at 5(0) and ending at e l27T ^ d S(0) defines an open 
path in SU (d) and a topologically nontrivial closed path 
R(t) £ Adj(d). If this cyclic evolution were composed d 
times, we would get a trivial path in Adj(d), so that the 
number of nonequivalent classes is given by d. 
The total phase can be written as 

<j> tot = axg{Tr[a^(0)a(t)]} = <\> A + <j> B 

+ axg{Tr[a^O)U A (t)a(0)Ul(t)]} , (3) 

while the dynamical phase is, 



^dyn 



dt' Tr[a f (t')a(t')} = 4> A + ( j) I 



i \ dt' Tr[p B {0)U A U A + P \{Q)UlUl 



(4) 



where p A — (S jf Q 2 S)*, p B — Q 2 - For cyclic evolutions we 
have U 3 {t) = e l2vrn ^ d U 3 (0). Then, the total generated 
phase is 4>tot — (pA + 4*b + 2irn/d, n = n A + 71b, where 
the values n ^ 0, d,2d, correspond to topologically 
nontrivial paths. As already discussed, the total phase is 
always written as a dynamical plus a geometric part. In 
order to consider a fractional phase as topological, it must 
be built only as a part of the geometric phase, receiving 



no relevant contribution from the dynamical part. This 
means that at any time t, < t < r, we must have, 



dt'Tr[p B (0) U\U A + p\(0) UlU* B ] = . (5) 



This is satisfied by the maximally entangled states, for 
every possible local evolution Uj. In this regard, the 
invariant quantities in the evolution can be written as 
Tr[(Q 2 ) p }. In terms of the concurrence we can write 



Q 2 = {l/d)I+^C 2 n -C 2 q.f 



(6) 



where C m — y / 2(d — l)/d. The C — value corresponds 
to separable states. For maximally entangled states 
C = C m , giving Q 2 = (l/d)I, and p A = Pb = {\/d)I. 
In addition, for any Uj £ SU(d), the matrices UjUj are 
combinations of the generators T a . Therefore, using this 
information, the trace in the integrand of eq. ([5]) van- 
ishes. 

Now, let us consider an evolution on the first qudit 
A. In this case, (ip(0)\tp{i)) = Tr[Q 2 (0) U A (t)], while the 
dynamical phase is, 



dt'Tr[Q 2 (0)U[(t')U A (t')} . (7) 



These phases do not depend on 5(0) so that for simplicity 
we can consider 5(0) = /, that is, U A (t) = S(t). For 
qubits T a = <j a /2 (a = 1,2,3), where a a are the Pauli 
matrices. We shall assume that the basis is chosen so 
that q(0) = e 3 , that is, Q 2 = 1/2 + y/l - C 2 cr 3 /2. The 
unitary sector of the state evolution can be put in terms 
of Euler angles, in the form U A (t) = U m (t) Vz{t), where 



U„ 



e -i V T 3e i6T 2e i v T 3 



V 3 



JxT 3 



(8) 



Note that for cyclic evolutions, U m (0) — U m (r), while 
V 3 (0) = ±V 3 (t). In addition, U m can be expanded in 
terms of /, T\ and T2, as the term proportional to T 3 is 
obtained from Tr[T 3 e~^ T3 e ieT2 e^ T3 ] = Tr[T 3 e i9T2 } = 
0. Here, we have used that the latter exponential is a 
combination of / and Ti. Using a similar expansion for 
e J ^ T 3, we arrive to the conclusion that the terms in U m 
proportional to Ti, T 2 do not contribute to (ip(0)\ip(t)}. 

With the ingredients above we can work out the ex- 
pression for the time evolving overlap 



(i>(Q)\ip{t)} = e i,pA cos- 



cos — + i\J 1 — C 2 sin — 
2 2 



In terms of Q 2 , U m , and V 3 , the dynamical phase is 



(9) 



•) A - 1 



rt 1 

i / dt' -Tr[(I+Vl-C 2 a 3 ) 



x {UlU m + VjV 3 )} 



(10) 
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FIG. 1. Time evolution of the quantum state overlap for a 
pair of qubits with different concurrences. 



Using V3 = i (x/2) 03 V3 and defining the unit vectors rh a 
so that U m &a Um = ™a, • cr , we get 



= arctan 



y/l - C 2 tan( X /2) - v 7 ! - C 2 ( X /2) 



+ Vl-C 2 (*/2), 



(11) 



with $ = J Q (it' mx • ^2 • In the last term, the frame rh a 
depends on G [0, 7r) and tp £ [0, 27r] defining a point 
on 5 2 , the surface of a sphere with unit radius. Then, 
iha(0,(p) is a mapping S* 2 — >• m a , and the evolution on 
this sector is given by a curve, defined by 0(t), (p(t), con- 
tained on S 2 . In this regard, for a cyclic evolution, one 
easily shows that $ = J7, where fi is the solid angle sub- 
tended by the closed path [15j, |l6j. This term can be 



associated to the usual Berry phase for a single qubit. 

For a general evolution, we see that for product states 
(C = 0), the first two terms in eq. (fTTj) cancel each other 
while the last term coincides with the one given by the 
usual picture of the Bloch sphere evolution of a single 
qubit. On the other hand, for maximally entangled states 
(C = 1), the last two terms vanish while the first term 
can assume only two discrete values or tt. In fig Q] this 
evolution is represented as paths in the complex plane, 
where the overlap (?p(0)\ip(t)) is plotted for different val- 
ues of the concurrence. This path degenerates to a circle 
for product states and to a straight line on the real axis 
as the concurrence approaches its maximum value (7=1. 
It gives a graphical picture of the phase jump between 
and 7r discussed in ref. [Ioj |. This jump occurs when the 
evolving state crosses the subspace orthogonal to the ini- 
tial one. Note that the solid lines in figQ] correspond to 
closed paths since points P and P' represent physically 
equivalent quantum states. Dashed lines correspond to 
additional closed paths. 

Now, let us study a simple nontrivial path that gener- 
alizes the l^-sector for qubits (cf. eq. ©) to the case of 



qudits. Consider an evolution of the form {/^(t) = Vjv(i), 
V N (t) = e^ E , X (0)=0, (12) 
where E is a diagonal traceless matrix with components, 



E — 



(1/d), a=l,...,d-l 
(1/d) - 1, a = d. 



(13) 



This matrix can be written in terms of the iV'th generator 
of SU(d), N = d 2 - 1: E = C m T N . When x (r) = 2tt, it 
is simple to see that Vm{t) 
q(0) = £n, we have, 



rJ2 



I. In the case where 



! (0) = (1/d) I + ^l-(C/C m ) 2 E 



(14) 



By expanding the exponential in eq. (|1 2|) and using 
eq.flUJl we get, 



(15) 



with A = ^± + -C 2 and B = 1 — A. Using 

Vfr Vn — i X E in the dynamical phase, we arrive at 



= arctan 



.4 sin ■ 



£?sin 



(i-rf)x 



„4cos ^ 



Bcos 



d 



Vc 2 



In the above example, for maximally entangled states, 
when d > 3 the total phase changes continuously from 
to 2n/d, and the evolving state never becomes orthogonal 
to the initial state. This is in contrast to what happens 
in the d = 2 case. The minimum value for | (i/j(0)\ip(t)) | 2 
is (A - B) 2 = (^) 2 , attained when x = T- For d = 3, 
the minimum overlap is (1/3) 2 . 

It is interesting to look for topologically nontrivial evo- 
lutions for qudits with similar properties to those dis- 
played by qubits. In the d — 3 case, this can be realized 
as follows. Let us consider the path J7a(x(*))> contin- 
uously evolving from Ua(0) = I to Ua(2tt) = e l2Tr / 3 /, 
defined by a diagonal unitary matrix with nontrivial ele- 
ments e** a such that fa = 2 X /3 + [2(tt - C)/3]©(x ~ O) 
fa = — 2 X /3, and fa = —{fa+fa); O(x) is the Heaviside 
function. For maximally entangled states, we have 



2cos(f )] , 
2cos(^±z0) 



e 3 



xe [0,tt] 

X e [tt, 2tt] . 



Then, we see that the total phase vanishes in the first 
part of the evolution, while it takes the fractional value 
27r/3 in the second part. In addition, at x = 7r, when the 
phase changes discontinuously, the state \ip(t)) becomes 
orthogonal to the initial state. 

Both qutrit evolutions are represented in figl2^, where 
the overlap (ip(0)\ip(t)) is plotted in the complex plane 
for maximal concurrence. The first cyclic evolution from 
P to P' is represented by the solid black line clearly show- 
ing that the overlap between the initial and the evolving 
quantum states never vanishes. On the other hand, the 
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FIG. 2. (a) Complex plane representation of the quantum 
state overlap for a pair of qutrits with maximal concurrence. 
Two different time evolutions are considered, (b) The corre- 
sponding stepwise evolution of the geometric phase. 



second evolution (red online) shows a path crossing the 
origin of the complex plane, where the evolving quantum 
state becomes orthogonal to the initial one. The dashed 
lines correspond to additional closed paths defining three 
vertices which evidence the fractional phase values. In 
figEb, we plot the associated geometric phase evolution, 
showing a stepwise behavior with two jumps between the 
fractional values 0, 2w/3, and Att/3. For the first evolu- 
tion (black) smooth jumps occur, while for the second 
evolution (red online) they are discontinuous. 

As a conclusion, in this letter we studied unitary local 
operations on a pair of qudits, showing that fractional 
phases naturally appear when cyclic evolutions are con- 
sidered. These fractional values are related to different 
homotopy classes of closed paths in the two-qudit Hilbert 
space. The geometric phase has been calculated in terms 
of the I-concurrence introduced in ref.[13|. In the case 
of maximally entangled states, the fractional values orig- 
inate solely from the geometric part of the phase evolu- 
tion, since the dynamical part vanishes at all times. 

The fractional phase of maximally entangled states is 



built in a stepwise evolution, where the phase jumps be- 
tween discrete values in steps of 2w/d. For qubits this 
jump is strictly discontinuous, while for qutrits, it may be 
discontinuous or not, depending on the particular evolu- 
tion considered. Due to its stepwise evolution, we expect 
the fractional phase acquired by maximally entangled qu- 
dits to be particularly robust against the influence of the 
environment. In order to produce a relevant change, any 
external noise would have to cause a large fluctuation, 
driving the two-qudit system through a phase step. Since 
the phase jump for qubits is strictly discontinuous, its ro- 
bustness should be even more pronounced. These results 
can be important to proposals of quantum gates based 
on topological phases. 
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